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ABSTRACT 



Neutral nickel (Ni I) is abundant in the solar atmosphere and is one of the 
important elements that contribute to the emission and absorption of radiation 
in the spectral range between 1900 and 3900 A. Previously, the Solar Radiation 
Physical Modeling (SRPM) models of the solar atmosphere considered only few 
levels of this species. Here we improve the Ni I atomic model by taking into 
account 61 levels and 490 spectral lines. We compute the populations of these 
levels in full NLTE using the SRPM code and compare the resulting emerging 
spectrum with observations. 

The present atomic model improves significantly the calculation of the solar 
spectral irradiance at near-UV wavelengths that are important for Earth atmo- 
spheric studies, and particularly for ozone chemistry. 

Subject headings: Radiative Transfer, Ni I Line Profiles, UV, Solar Spectral Irradiance 
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Introduction 



The radiative flux spectrum incident at the top of the Earth atmosphere is known 
as the solar spectral irradiance (SSI). The near-UV range of the SSI plays a fundamental 
role in the ozone chemistry. Most of the Earth atmospheric ozone is produced in the 
stratosphere as a result of photo-dissociation of molecular oxygen by UV Solar radiation in 
the Schuman-Runge bands at wavelength shorter than ~ 1750 A, a nd destroyed by ozon e 



photo dissociation in the Hartley-Huggins bands (1900 to 3900 A) (IFontenla et al.lll999 ). 



A change in the amount of near-UV SSI affects directly the destruction of ozone leading to 
tropospheric effects. Ozone changes al so affect the am ount of biologically harmful near-UV 



radiation reaching the Earth's surface (IBrasseur 



19971 ) . For these reasons, there is a need for 



an accurate calculation of the solar radiation that reaches the top of the Earth atmosphere 
in the near-UV range. 

The Solar Radiation Physical Modeling (SRPM) system is a code library that permits 
one to compute the radiance (emergent intensity from the solar surface) spectrum of the 
Sun based on physical models of the solar atmosphere (see Fontenla et al. 2006, 2007, and 
2009). The code solves the radiative transfer and non-LTE statistical equilibrium assuming 
a specified one dimensional plane— parallel geome try or a spherical symmetry depending 



of the study. The paper by 



Fontenla et al. 



( 1201 ll ) applied the SRPM code and a set of 



physical models of various features observed on the solar surface to the calculation of high 
spectral-resolution SSI from extreme ultraviolet to far infrared. The results obtained in 
that paper allow one to study the effects of SSI changes on the Earth and other planetary 
atmospheres. Comparisons show a generally very good agreement between the observed 
and synthetic irradiance, but significant differences occur at some wavelengths. 

Important discrepancies, especially in the near-UV range, reflect the need for more 
accurate and complete atomic models of some important species emitting and absorbing in 
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that range. Ni I is an important species, considering its abundance in the solar atmosphere, 
and has strong lines in the near-UV range. The present paper focuses on improving 
the calculation of the solar spectra in this range by improving the Ni I atomic model 
and using the SRPM system to compute the non-LTE populations and the emergent 
intensities. Comparisons are shown between the new results and observations of radiance 
and irradiance. 

Only a few publications concerning Ni I NLTE calculations are found in the literature, 
most of them about the formation of the Ni I 6769.64 A line. In helioseismology, this line 
is used for observing line-of-sight velocity (Doppler shift) by the Michelson Doppler Imager 
on-board SOHO and by the ground-based Global Oscillations Network Group. 



Brulsl (11993h studied the Ni I 6769.64 A line using an ato mic model with 



levels. The atomic data used for the calculation were that by 



Corliss fe Sugar 



9 energy 



(1981) 



The level populations and emer gent intensities were computed in NLTE using the VAL C 



model by 



Vernazza et al. 



( 11981I ). This paper confirmed that the line is a good choice for 



helioseismology. 



More recently, 



Norton et al. 



(12006f ) compared the Ni I 6769.64 A and the Fe I 6173 
A lines, both candidates to be used by the Helioseismic and Magnetic Imager (HMI) 
on-board the Solar Dynamics Observatory (SDO) for observing Doppler velocities and 
magnetic fields. They adopted a Ni I atomic model with 25 energy levels and the same 
source of atomic data as in the pre viously mentioned wor k. The NLTE calculations were 



carried out again using model C by 



Vernazza et al. 



model M by 



Maltbv et al. 



(1198ll ) and also using sunspot umbra 



( 119861 ). The HMI team selected the Fe I spectral line because of 



its better performance for magnetic diagnostics. 



Considering the spectral lines and term system of Ni I, the results of 



( 119931 ) increased the number of known Ni I levels reported by 



Litzen et al. 



Corliss fc Sugar 



( Il98lh from 
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185 to 286, and the number of classified Ni I lines were expanded from 1071 to 1996. This 
work is the source used by the National In stitute of Standards and Technology (NIST) 



database version 4 (IRalchenko et al 



20111 ). The lowest 61 energy levels in the NIST 
tabulation are adopted in the present paper because they correspond to the most important 
lines particularly in the near-UV. 



In the present work we calculate the 



solar features described by 



Fontenla et al. 



Quiet Sun spectrum, as defined by the mix of 



(1201 if ), using the same set of physical models but 
improving the Ni I atomic model from the original 10 to the present 61 levels, and compare 
the computed spectrum with observations of several Ni I line profiles. 

This paper is arranged as follows: the new Ni I atomic model is presented in Section 2. 
In Section 3 we show some results concerning the radiative transfer and NLTE calculations. 
In Section 4 we compare the computed quiet sun spectra with observations. Finally, in 
Section 5 we discuss our results. 



2. The improved Ni I atomic model 

The improvement of the Ni I atomic model is primarily motivated by producing a 
better computation and comparison with observations of the near-UV SSI. 

Observations in this spectral range are scarce because the atmospheric absorption, 
especially by ozone, make ground observations impossible. SSI observations from space 
usually have limited spectral resolution. 

A limitation for computing an accurate spectrum is the availa bility of accurate atom ic 
data for NLTE calculations. The NIST atomic spectra database, ( IRalchenko et al.ll201ll ). 
lists 503 Ni I spectral lines in the range 1900 to 7000 A, with 266 belonging to the near-UV 
range (1900 to 3900 A). 
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To reproduce most of the the strong Ni I lines in the near-UV rang e, it is necessary 



to con sider an atomic model with 61 energy levels, and as described by 



Ralchenko et al 



( 120 111 ) the corresponding sublevels defined by the J quantum number. Using this number 



of levels only 34 lines of the NIST listed lines are not included. 

The level energies, line wavelengths, and weighted oscillator strength gf values were 
taken from the NI ST atomic spe ctra database. The collisional strength parameters were 



obtained from the 



Seatonl (119621 ) approximation, and the photoionization cross sections 



were estimated from similar atomic structure species. The reason for these approximations 
is the absence of better Ni I atomic data in the literature. For the r adiative, Stark, and Van 



der Waals broadening parameters, we used the values estimated by 



Kurucz 



jl994h . 



In Table 1 we present the description of the 61 atomic levels and sublevels of energy 
considered. The first column lists the configuration, and the second column the associated 
term. The third column shows the values of the J quantum number that define the 
sublevels, and the fourth column presents the sublevel energies. The fifth column lists the 
level number in our scheme. 



3. The NLTE calculations 



We carried out the c alculations using the latest semi-empirical models of the solar 



atmosphere described by 



Fontenla et al 



(1201 ll ). The Quiet Sun spectrum was built by 



weighted averaging of the three models listed in Table 2 and shown in Figure 1. To 
reproduce the low solar activity state during the period 2008-2009, we considered the solar 
disk as covered by 80% of internetwork feature (model 1001, B), 19% of network feature 
(model 1002, D) and 1% of active network (model 1003, F). 



The populations for the atomic species listed in Table 2 of 



Fontenla et al. 



( 1201 lh were 
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all calculated in full NLTE as in iFontenla et al.l (1201 if ), but now with the new Ni I atomic 



model that increased the number of energy levels from 10 to 61. 

We calculated the LTE departure coefficients {b values) defined by the ratio of the 
computed level populations relative to those resulting from the Saha-Boltzman distribution 
to the next ion stage (Menzel definition, Rutten 2003). Figure 2 shows the calculated 6's 
for many levels that resulted from the atmospheric model 1001. 



Fontenla et al 



( 120091 1 . At heights 



The general behavior is similar to that described in 
above ~ 900 km (ie, above the temperature minimum region) the b values are greater than 
unity, with the largest value corresponding to the ground level. The opposite occurs in the 
low chromosphere and the temperature minimum region where the b values are smaller 
than unity. In general, as the energy of the level increases the departure coefficients are 
closer to unity, but fluctuations occur depending on how the various levels are connected by 
permitted transitions. 

In Figure 3 we plot the b\ departure coefficients of the ground level for the three 
models used to reproduce the quiet sun spectrum. Although the general behavior is similar, 
the region where b\ is smaller than unity is thinner, and closer to unity, as the minimum 
temperature of the atmosphere increases. 

In Figure 4 we show the ion density fractions of Ni I/Ni tota ;, Ni II/Ni to4a j, and Ni 
III/Ni tota z with height. Throughout the atmosphere the total Ni II density is nearly the 
total Ni density. Above ~300 km Ni I/Ni tota ; drops slowly. Above ~600 km Ni III/Ni tota ; 
begins to increase, although below the chromosphere-corona transition region it remains 
several orders of magnitude lower than the other two density fractions. 

In Figure 5, we plot the line source function (solid line) and the local Planck function 
(doted line) versus the optical depth at line center r for four important line transitions. 
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In a first approximation, line center forms at r = 1 for disk center and r = 2/3 for the 
irradiance. Considering this plot, departures from LTE produce a large difference between 
these functions for most of the plotted r values , and specially around r = 1. The heights 
at which r = 1 are as follows; ~ 597 km for the 2747.554 A line, ~ 945 km for the 3415.74 
A line, ~ 840 km for the 3567.39 A line, and ~ 323 km for the 6769.64 A line. Thus, the 
centers of these lines form in the low chromospheric layers at and below the temperature 
region where NLTE effects are strong, and their source functions cannot be well described 
by any simple approximation. 



4. Comparison of the calculated spectra with observations 



We carried out the calculations of the solar spectrum with the new Ni I atomic model. 



using the latest semi-empirical models of the solar atmosphere described by iFontenla et al 



(120111 ). The Quiet Sun spectrum was constructed by weighted averaging of the spectra from 
three models listed in Table 2 and shown in Figure 1. To reproduce the low solar activity 
state during the period 2008-2009, the solar disk is considered covered by the relative areas 
of features indicated in Table 2. 

The information about the Ni I lines and energy levels that we compare with 
observations are presented in Table 3. The last column in this Table indicates the figure 
number where this transition is plotted. Many of the line profiles shown in the Figures 
correspond to transitions between energy levels that were never before computed in full 
NLTE. 

Computed spectral lines in the range from 3290 to 7000 A are compared with the 



observations at disk-center in the FTS Solar Atlas by 



Neckell (119991 ). The absolute intensities 



of these observations were from the Kitt Peak FTS-data as described by Neckel and Labs 
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1984a, b) and Neckel (1994) . The absolut e calibration of this disk-center spectrum is not 



very reliable as explained by 



Neckell (ll999[ ). For this reason, we scaled the observations to 



match the continuum near the line profiles with the calculations. Figures 6 to 8 show the 
comparison of several line profiles and the scaling factors. 

The last panel in Figure 7 shows the Ni I 6769.64 A (vaccum wavelength) line used in 
SOHO/MDI and GONG instruments. 

The differences between the two computed profiles plotted in these Figures are due 
to the different number of levels in both NLTE calculations that was mentioned before. 
Several Ni I lines were absent in the spectra calculated by Fontenla et al. (2011), while 
others are too deep (e.g. dashed trace in Figure 7). The later are improved in our present 
calculation by the inclusion of many more levels. 

In some cases the calculated line profile is wider than the observed one. We also 
recomputed the line profiles usi ng the broadening parameters included in the Vienna 



Atomic Line Database (VALD) (IKupka et al. 



20001 ) . but the results obtained were similar. 



Comparisons with observations were also carried out with near-UV irradiance spectrum 
in the range from 190 to 310 nm. The data from space instruments by Thuillier et al. 
(2003), and by Harder et al. (2010) have a reliable absolute calibration. However, these data 
have low spectral resolution. In Figure 9 we show the compari son between our comp uted 



Thuillier et al. 



(120031 ). To 



spectral irradiance and the Composite 3 reference spectrum by 
simulate the resolution of the observatio ns, the calculate d spec tral irradiance was convolved 



with a 1 nm F WHM cos 2 function as in iFontenla et al. 



(1201 lh . 



The observa t ions o f the irradiance spectrum by 



Hall fc Anderson! (119911 ) and 



Anderson fc Halll ( 119891 ). hereafter H&A, were obtained from stratospheric balloon flights 



that are very difficult to calibrate because of the atmospheric absorption. However, these 
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observations have a fairly good spectral resolution of 0.01 nm. For the present comparison 
the H&A spectrum was scaled to overall match the values in the calculated spectrum. In 
Figure 10 we compare several computed Ni I line profiles with the scaled H&A irradiance 
observations. The comparison shows good agreement between calculations and observations. 



5. Summary 



This work improves the Ni I atomic model to include the majority of the lines in the 
near UV range that are listed in the NIST atomic database. Several of these lines are very 
strong and have not been calculated in full NLTE by previous papers. 

The results shown here display a fairly good agreement with the observed Ni I lines 



because we include full NLTE calculations for an atomic model with 61 leve 
fine st ructure sublevels. This agreement does not require any changes of the 



s, which have 



Fontenla et al 



( I201ll ) atmospheric models. 



Also, we reproduce well the Ni I 6769.64 A used in MDI and GONG observations, 
allowing a better understanding of the line formation in different solar features, and 
investigating the impact of solar activity on its behavior. 

The results presented in this work stress the importance of full NLTE calculation and 
the use of a sufficient number o f levels for reali s tic co mputation of line profiles. From these 
results it is clear that using the 



Fontenla et al. 



( 1201 ll ) models for spectra calculations with 



LTE codes will not properly compute the deep Ni I lines and likely others. This is because 
the calculations here show that the line source function substantially departs from LTE in 



the low chromosphere, and at t 
produce wrong results with the 



re temperature minim um. LTE values will systematically 



Fontenla et al 



(1201 lh models. 
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Table 1: The first 61 levels of Ni I used in the present paper (IRalchenko et al. 



Configuration 


Term 


J 


Level (cm x ) 


Level number 


3d 8 ( 3 F)4s 2 


3 F 


4 


0.000 


1 






3 


1 332.164 








2 


2 216.550 




3d 9 ( 2 D)4s 


3 D 


3 


204.787 


2 






2 


879.816 








1 


1 713.087 




3d 9 ( 2 D)4s 


X D 


2 


3 409.937 


3 


3d 8 ( 1 D)4s 2 
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4 


3dl0 


l S 
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5 


3d 8 ( 3 P)4s 2 
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6 
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2 
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6 


27 260.894 


9 






5 


27 580.391 








4 


28 068.065 








3 


28 578.018 








2 


29 013.206 





20111 ). 
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Table 1: The first 61 levels of Ni I used in the present paper (iRalchenko et al 



20111). 



Configuration 


Term 


J 


Level (cm l ) 


Level number 


3d 8 ( 3 F)4s4p( 3 P°) 


5po 


5 
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10 






4 


29 084.456 
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1 
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3d 9 ( 2 D)4p 


3po 
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1 
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12 






3 
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13 






2 
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14 
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15 






4 
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3 


31 786.162 




3d 9 ( 2 D)4p 


lpo 


3 


31 031.020 


16 
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Table 1: Continuation 



Configuration 


Term 


J 


Level (cm x ) 


Level number 
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Table 1: Continuation 
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Table 1: Continuation 
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Table 2: Quiet Sun feature components and their respective models. We follow the feature 



and model index designations as in Fontenla et al. (2011). 



Feature 


Description 


Photosphere-Chromosphere 
Model Index 


Relative area on 
the Solar Disk 


B 


Quiet Sun inter-network 


1001 


80 % 


D 


Quiet Sun network lane 


1002 


19 % 


F 


Enhanced network 


1003 


1 % 
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Table 3: Transition lines plotted in Figures 6, 7, 8 and 10. In column 1, the vacuum wave- 
length. In column 2 and 3, the energy value for the upper and lower sublevel respectively 



involved in the transition (IRalchenko et al. 



20111). 





Upper level (cm x ) 


Lower level (cm x ) 


Figure number 


2296.69 


43654.974 





10 


2326.515 


44314.904 


1332.164 


10 


2747.554 


36600.791 


204.787 


10 


2822.122 


35639.122 


204.787 


10 


2982.516 


34408.555 


879.816 


10 


3003.36 


33500.822 


204.787 


10 


3370.531 


29668.918 





6 


3393.957 


29668.918 


204.787 


6 


3415.744 


29480.989 


204.787 


6 


3434.54 


29320.762 


204.787 


6 


3447.247 


29888.477 


879.816 


6 


3459.451 


30619.414 


1713.087 


6 


3493.956 


29500.674 


879.816 


7 


3516.057 


29320.762 


879.816 


7 


3525.544 


28569.203 


204.787 


7 


3567.39 


31441.635 


3409.937 


7 


3859.391 


29320.762 


3409.937 


7 


4717.081 


49777.569 


28578.018 


7 


4830.372 


49271.540 


28569.203 


8 


4919.737 


51306.038 


30979.749 


8 


5036.766 


49174.770 


29320.762 


8 


5081.949 


49158.480 


29480.989 


8 


5478.426 


32982.260 


14728.84 


8 


6769.64 


29500.674 


14728.84 


8 
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Fig. 1. — Atmospheric models for the solar features used to calculate the quiet sun spect 
These models were built by Fontenla et al. (2011) 
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Fig. 2. — Departure coefficients b for selected levels of Ni I in the atmospheric model 1001. 
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Fig. 3. — Departure coefficients bl for the ground state of Ni I in the atmospheric models 
used to calculate the solar spectrum 
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Fig. 4. — Ni I, Ni II and Ni III densities relative to Ni totai density for the atmospheric model 
1001 as a function of height. 
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Fig. 5. — Line Source Function (S, solid line) and Planck Function (B, doted line) as a 
function of the optical depth at line center, r, for selected spectral lines. 
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Fig. 6. — Comparison between observed and synthetic line profiles at disk center. Solid lines 
show our calculation, dashed li nes the spect ra calculated by Fontenla et al. (2011), and 



doted lines the observations by 



Neckell (Il999f ). The scale factors, in increasing wavelength 



order, are 1, 1, 0.91, 1, 1.18, 1.18. All wavelengths are in vacuum. 
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Fig. 7. — Comparison between observed and synthetic line profiles at disk center. Solid lines 
show our calc ulation, dashe d lines the spectra calculated by Fontenla et al. (2011), and 



doted lines the 



Neckell (119991 ) Solar Atlas. The scale factors, in increasing wavelength order, 



are 1.28, 1, 0.91, 1, 1.18, 1.25. All wavelengths are in vacuum. 
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Fig. 8. — Comparison between observed and synthetic line profiles at disk center. Solid 
lines show our calc ulation, dashe d lines the spectra calculated by Fontenla et al. (2011), 



and doted lines the 



Neckell (119991 ) Solar Atlas. The scale factor for all the lines is 1.28. All 



wavelengths are in vacuum. 
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Fig. 9. — Comparison between observed and synthetic spectral irradiance between 190 and 
310 nm. S olid lines show o ur ca lculation, dashed lines the Composite 3 reference spectrum 
derived by iThuillier et all (hooa ) 
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Fig. 10. — Comparison between observed and synthetic irradiance line profiles. Solid lines 
show our calculation, and doted lines the H&A observations. All wavelengths are in vac- 
uum. The scale factors between the observed and the calculated profile are, in increasing 
wavelength order, 1.9, 3.3, 3.5, 1.8, 0.75 and 0.83 



